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ABSTRACT 

White dwarfs inspiraling into black holes of mass Mbh ^ lO^M© are detectable sources of 
gravitational waves in the LISA band. In many of these events, the white dwarf begins to lose 
mass during the main observational phase of the inspiral. The mass loss starts gently and can 
last for thousands of orbits. The white dwarf matter overflows the Roche lobe through the Li 
point at each pericenter passage and the mass loss repeats periodically. The process occurs 
very close to the black hole and the released gas can accrete, creating a bright source of radi- 
ation with luminosity close to the Eddington limit, L ~ 10^^ erg s~^. This class of inspirals 
offers a promising scenario for dual detections of gravitational waves and electromagnetic 
radiation. 
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1 INTRODUCTION 

One of the goals of the Laser Interferometer Space Antenna (LISA) 
mission is to detect gravitational waves from compact stellar ob- 
jects spiraling into massive black holes, a class of events called 
extreme-mass-ratio inspirals (EMRIs) (e.g., Hils & Bender 1995; 
! Gair et al. 2004; Barack & Cufler 2004; see Hughes 2009 and 
Sathyaprakash & Schutz 2009 for recent reviews). Of particular in- 
terest are sources of coincident gravitational and electromagnetic 
, radiation. Besides providing unique information on the nature of 
' the event, such dual detections will lead to a new version of Hubble 
diagram that is based on the gravitational distance measurements 
(e.g. Bloom et al. 2009; Phinney 2009). 

Inspirals into black holes of masses Mbh ^ (10^ — 1O^)M0 
produce gravitational waves in the frequency band where LISA is 
most sensitive. Normal stars are tidally disrupted well before they 
approach such black holes, and therefore discounted as possible 
LISA sources (however, see Freitag 2003). Inspirals of compact ob- 
jects are guaranteed sources of gravitational waves, however most 
of them are not promising for dual detections. In particular, stellar- 
mass black hole or neutron star inspirals are not expected to gen- 
erate bright electromagnetic signals. Only inspiraling white dwarfs 
(WDs) offer a possibility for dual detection (Menou et. al. 2008; 
Sesana et al. 2008). WDs can be tidally disrupted very close to the 
black hole and then could create a transient accretion disc with Ed- 
dington luminosity. 

Estimated EMRI rates are high enough for observations with 
LISA (e.g. Phinney 2009). The expected fraction of WD inpirals 
among all EMRIs depends on the degree of mass segregation in 
galactic nuclei, which favors stellar-mass black holes over white 
dwarfs in the central cluster. The abundance of white dwarfs also 



depends on the details of stellar evolution, which are not completely 
understood. A fraction of WD inpirals as large as ^ 10% has been 
suggested (e.g. Hopman & Alexander 2006a). 

The orbital parameters of WD inspirals are also uncertain. 
EMRIs form when a compact object is captured onto a tight orbit 
whose evolution is controlled by gravitational radiation rather than 
random interactions with other stars in the central cluster around 
the massive black hole. Two main channels exist for EMRI for- 
mation: capture of single stars and capture of binary systems (e.g., 
Hils & Bender 1995; Sigurdsson & Rees 1997; Ivanov 2002; Gair 
et al. 2004; Hopman & Alexander 2005, 2006a,b, 2007; Miller et al. 
2005; Hopman 2009). In the single-capture scenario, the shrinking 
orbit can retain a significant eccentricity until the end of inspiral. In 
contrast, the binary-capture scenario leads to nearly circular orbits 
(Miller et al. 2005). 

In this paper, we focus on WD inspirals that are not completed 
because the star is tidally disrupted before its orbit becomes unsta- 
ble. We argue that the WD begins to lose mass very gently and, 
for thousands of orbital periods, this process resembles accretion 
through the Li point in a binary system rather than a catastrophic 
disruption. This offers a possibility of simultaneous observation of 
the inspiral by LISA and traditional, optical and X-ray telescopes. 
Previous work on tidal deformation of a WD orbiting a massive 
black hole focused on two extreme regimes: (i) weak deformation 
was studied analytically using perturbation theory (e.g., Rathore et 
al. 2005; Ivanov & Papaloizou 2007), and (ii) strong deformation 
leading to immediate disruption was simulated numerically (e.g. 
Kobayashi et al. 2004; Rosswog, Ramirez-Ruiz & Hix 2009). The 
regime considered in the present paper is different from both cases 
explored previously. It involves an extended phase of strong defor- 
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mation with small mass loss, which we call 'tidal stripping' below. 
The mass of the WD remains almost unchanged during this phase 
and continues to emit gravitational waves. Even a small orbital ec- 
centricity, e.g. e = 0.01, implies that tidal stripping occurs only 
near the pericenter of the orbit, during a small fraction of each or- 
bital period. Thus, the mass loss is expected to h^ periodic, possibly 
leading to a periodic electromagnetic signal from the inspiral. 



2 TIDAL STRIPPING 

2.1 Onset of mass loss 

Consider a WD orbit with semi-major axis a and eccentricity e. 
The orbital parameters gradually evolve as a result of gravitational 
radiation (Peters 1964) 
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where dot denotes the time derivative of the secular evolution (av- 
eraged over the orbit). Both a and e slowly decrease with time. The 
fractional change of the pericenter radius Vp = a(l — e) in one 
orbital period P = 2n{a^ /GMmf^'^ is 
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For the typical parameters of the problem considered in this pa- 
per, a ^ 10~^. Two inaccuracies in the above formulas should 
be noted: (i) Equations ([T])-(|3]) are valid only if Vp ^ GMbu/c^- 
They become approximate during the most interesting phase of the 
inspiral when rp is only a few rg = 2G'Mbh/c^. (ii) The equations 
neglect the effects of mass loss on the evolution of the orbitjj 

The mass loss begins when the tidal acceleration created by 
the black hole at the WD surface, (GMBH/rp)R, becomes com- 
parable to GM/R^, where R is the radius of the WD. A similar 
condition for the onset of mass transfer in synchronous binary sys- 
tems is that the donor star fills its Roche lobe. In the more com- 
plicated case of eccentric non- synchronous binary systems, mass 
transfer may be approximately described using an instantaneous 
effective Roche lobe (Sepinsky et. al. 2007), neglecting all the ef- 
fects beyond the quasi-static limit (e.g. Ritter 1988). We will use 
the following approximate condition for the mass-loss onset, 

where R is the radius of the unperturbed star, before it experiences 
any tidal deformation, and 7 is a numerical constant. The exact 
value of 7 depends on the mass ratio M/Mbh and the orbital pa- 
rameters a and e. Even more importantly, it also depends on the ro- 
tation of the WD and the history of its tidal heating, none of which 
is known. If Rt is interpreted as the effective Roche-lobe studied 
by Sepinsky et al. (2007), their results give 0.39 ^ 7 ^ 0.59, with 
7 = 0.49 for synchronous binary systems. Our conclusions are in- 
dependent of the exact value for 7. In numerical examples where it 
needs to be specified, we assume 7 ^ 0.5. 

^ Let 5M be the mass lost in one orbit. The maximum correction to a is 
~ 5M/M, and it is significant only when 5M exceeds aM. Section 2.3 
suggests that the value of a is unimportant at this stage, as the evolution 
becomes controlled by the mass loss itself, not the drift of the pericenter. 
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Figure 1. Radius ro where tidal stripping begins, in units of the 

Schwarzschild radius r^ = 2GMbh/c^ - 3x1O1O(Mbh/1O^M0) cm, 
is shown as a function of the WD mass. The three curves correspond to 
^bh/10^^0 = 0-5, 1 and 5. The value of ro was estimated using the 
simplified equation (|6), which is non-relativistic and neglects the effect of 
the black-hole spin as on the tidal force. The shaded region shows radii 

min depends on as and 
3rp for {as = 0, e = 0}, 
= rg/2 for {as = l,e = I}. 



r < rn 



where no stable bound orbits exist, 



orbital eccentricity e, in particular rmin 
^min = 2rg for {as = 0, e = 1} and r^: 



The orbit slowly shrinks due to gravitational radiation and 
condition ^ will be first met when the pericenter radius rp reaches 
the value ro estimated below. We adopt a simple model for the un- 
perturbed WD star: a non-rotating cold sphere of uniform chemical 
composition, supported by the pressure of degenerate electrons. It 
is straightforward to obtain numerically the mass-radius relation 
for such a star. We find that it is well approximated by the follow- 
ing formula (with less than 2 per cent error for O.2M0 < M < 
I.4M0), 
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where Mch = 1.43M0 is the Chandrasekhar mass, /3 = 0.447 and 
Ri, = O.O13i?0. Substituting R{M) in equation Q, one obtains 
the pericenter radius at which mass loss begins 
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Figure[T]shows ro/rg as a function of WD mass for different Mbh- 

2.2 Estimate for mass loss in one orbital period 

As the star begins to lose mass, it continues to do so repeatedly: 
at each pericenter passage the star overflows the Roche lobe for a 
short time. The hydrodynamics of the surface layers of the tidally 
stripped star is complicated. It is clear however that the mass lost 
in one pericenter passage, 6M, will be tiny for the initial mass loss 
episodes when the star surface barely touches the Roche lobe. The 
natural small parameter in the problem is A/R = (R — Rt)/R. 
The onset of mass loss is defined by A = and we assume that the 
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dependence of 6M on A/R <^ 1 can be expanded in series whose 
leading term has the form, 



SM _ fA 

W-^ [r 



A = R-Rt<^R, 



(7) 



where A is a numerical factor. 

For illustration, consider a toy model. Suppose that after the 
pericenter passage the star loses the surface shell A <^ R whose 
mass is estimated using the structure of the unperturbed star. 



SM = 47vR 
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Here z is the depth measured inward from the surface, and p{z) is 
the mass density of the unperturbed star at depth z. For simplic- 
ity, let us assume the poly tropic equation of state in the surface 
layer P = Kp^ with 7 = 5/3, the same as in the deeper region 
where electrons become degenerate. P and p change continuously 
between the degenerate and non-degenerate regions, and therefore 
K in the surface layers must be the same as for non-relativistic 
degenerate electron gas, 

1 /sy/^ K^ 

where /Xe ~ 2 is the mean molecular weight per electron. The hy- 
drostatic balance dP/dz = p GMjR^ gives 
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and evaluating the integral in equation ^ one finds 6M/M = 
A(A/i?)^ with 
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In this toy model, 6M = (6/5)[p{A)/p]{A/R) where p = 
3M/47vR^. The small factor (A/i?)^ < 1 in ^M results from 
the small thickness of the surface layer A and the small density of 
this layer, p{A)/p - (A/i?)^/^ 

Additional factors may enter a more realistic model. The du- 
ration of the mass-loss episode St may be so short that only a frac- 
tion of the A-layer is lost. 6t likely scales with some power of 
A/R. One could formulate a time-dependent model for the mass- 
loss episode by evaluating Rt along the orbit around the pericenter 
(instead of just one point r^). If one assumes that the mass-loss 
episode occurs where Rt < R, then St ^-^ {A/RY^'^t where r 
is the sound-crossing time of the star. A complete hydrodynamical 
model is three-dimensional as the mass loss is asymmetric: gas will 
flow through the Li point, as happens in close binary systems. The 
flow velocity may be comparable to the sound speed in the outer 
layers of the star (which is smaller than the sound speed in its in- 
terior). Perhaps future hydrodynamical calculations will give the 
duration of the mass loss episode and the resulting ^, A and 5M. 
We do not know the exact values of V^ and A and keep them as pa- 
rameters. The illustrative numerical example shown below assumes 
V^ and A given in equation ([TT1) . 

From equations ^ and ^ one finds 

A ^ ^ ^ ^ rp (mV^^ / Mch-Mo y 

R R ro \MoJ \Mci.-M ) ' ^ ^ 

For the first mass loss episode, M — Mo and r^/ro ^ 1 — a, 
which implies A/R r^ a ^ 10~^. Equation i\2\ shows that A 
grows with each passage of the pericenter as rp decreases (due to 
gravitational radiation) and M decreases (due to mass loss). 



2.3 Evolution of mass loss over many orbits 

N orbits after the star reached ro, the pericenter radius is given by 
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ro 
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aN <^ 1, 
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where a is given by equation ([3]). Here we assumed that rp de- 
creases by arp in each orbit, neglecting the effect of mass loss on 
the orbit (which is likely to be a poor approximation at late stages 
of mass loss, e.g. Bildsten & Cutler 1992). 

Let X be the mass fraction of the star that has been lost over 
N orbits. 
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As long as X <^ 1, one can expand A in x and aN and keep only 
the leading linear terms. 
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The term aN describes the change in Rt/Ro due to the decreasing 
rp while the small changes in M and R are neglected. The term 
Bx describes the effect of decreasing M on Rt and R while the 
small drift of the pericenter is neglected. 

It is convenient to treat N ^ 1 as a continuous variable 
and describe the mass loss by the differential equation dM/dN = 
—SM. Then substitution of equation jTSl to equation ([7]) gives the 
differential equation for x{N), 
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dN 
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One can see from this equation that there are two stages of mass 
loss: (a) Bx <^ aN and (b) Bx ^ aN. Assuming ip > 1, we find 
the solutions for the two regimes. 
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Equation JTsl l applies only as long as x <^ 1, however it gives an 
estimate for the number of orbits to complete disruption N^ . The 
values of A^i and N^ can be evaluated by matching x and dx/dN 
for the two solutions at A^ = A^i , 



N, = {t^y^\('-^y^ 
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At A^ < A^i the growth of x is caused by the decrease in rp at 
practically unchanged mass and radius of the star. After A^i orbits, 
the growth of x accelerates as it is now controlled by the decreasing 
M (and increasing R) while the change in rp has a negligible effect. 
Equations Ct]) and dJ) imply that SM (x A^"^ for N < Ni and 
SM (X {N^ - Ar)^/(i-^) for N > Ni. 

For example, consider the toy model described by equa- 



tion (fTTT) . In this case, Ni ^ a' 



-3/5 



lO'' and N^ 



Ni = 



(4/21)A^i. The mass fraction lost after A^i orbits is x(A^i 



aNi/B 



v2/5 



10 . The detailed behavior of x{N) in the 



toy model is shown in Figure 2 for a WD with initial mass Mo 
O.GMq. The figure shows the solution of equation dx/dN 
A{A/R)'^ with A given by equation ([12]). 
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Figure 2. Lost mass fraction x = (Mo — M)/Mq after A^ orbits since the 
onset of tidal stripping. The WD has initial mass Mq = O.QMq. Different 
curves correspond to different eccentricities of the orbit; Mbh = 10^ M© 
is assumed for all cases. The insert shows log x plotted against log A^. 



2.4 Periodic mass-loss rate 

The mass-loss rate dM/dt is nearly periodic with the orbital period 
P. It is zero throughout most of the orbit and has a strong peak 
near the pericenter. To illustrate this behavior, we calculated the 
following, greatly simplified model. 



dM 
"dT 



{0- 



6M/r{M) 



Rt<R 
Rt> R 



(21) 



where R{M) is the radius of the unperturbed star of mass M, Rt 
is calculated everywhere along the orbit according to equation (|4]). 
SM{t) is given by equation (|7]) (it is evaluated using the local value 
of Rt); r — {Gp)~^^'^ is the sound-crossing time-scale of the star. 
Figure [3] shows the numerical solution of equation (|2TI ) for the 
last 35 orbits before disruption. The WD is assumed to have an 
initial mass Mq — O.6M0 and orbital eccentricity e = 0.9; the 
orbital period is P ?^ 2 x 10^ s. We plot the instantaneous mass- 
loss rate \dM/dt\ versus time tioss shown by a clock that ticks only 
when dM/dt / 0. The periodic peaks in \dM/dt\ coincide with 
the pericenter passages. 



3 DISCUSSION: ELECTROMAGNETIC COUNTERPART 

The standard picture of EMRI envisions an orbit that gradually 
shrinks due to gravitational radiation until its pericenter rp reaches 
rmin where the orbit becomes unstable and plunges into the black 
hole. The main observational phase of the inspiral is when rp de- 
creases from ^ 2rmin to rmin. For WDs with mass M ^ Mq 
inspiraling into a black hole with Mbh ~ 10^ M© the main obser- 
vational phase lasts ^ 10^ orbital periods, which may be compa- 
rable to one year, depending on the orbital eccentricity. We argued 
in Section 2 that the inspiraling WDs can experience an extended 
period of slow mass loss during observations by LISA. 

To summarize, the tidal stripping begins very gently because 
the pericenter of the orbit drifts inward slowly, by a tiny fraction 
a ^ 10~^ in one orbital period. This leads to many repeated 
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Figure 3. Instantaneous mass-loss rate M = \dM/dt\ during the last 35 
orbits before disruption. The WD with initial mass Mq = O.GMq is orbit- 
ing a black hole with Mbh = 10^ Mq; the orbit has eccentricity e = 0.9. 
The horizontal axis shows 'time during mass transfer' , which increases only 
when M / 0. The area under each peak is the mass lost per pericenter 
passage; it follows a power law with the number of orbits left to complete 
disruption, SM oc (A^^ — A^)~^-^, close to the results of Section 2.3. 



episodes of small mass loss. The process occurs at radius ro (eq.[6l), 
comparable to rmin- The star is strongly deformed by the tidal 
forces near ro but it barely touches its Roche lobe for a small frac- 
tion of the orbital period. As a result, the star loses a small amount 
of mass through the Li point at each pericenter passage. To our 
knowledge, this regime of tidal stripping was not explored by di- 
rect hydrodynamical simulations. Our crude estimates suggest two 
phases of stripping: the first phase lasts A^i ^ 10^ orbits until 
the star loses ^ 1 per cent of its mass. Then the mass loss ac- 
celerates: the decrease in M (and the corresponding increase in 
the WD radius) implies that the star overfills the Roche lobe more 
and more with every pericenter passage. As a result, the star loses 
the remaining 99 per cent of its mass in a few hundreds of addi- 
tional orbits. The mass lost in each individual episode during these 
last hundreds of orbits behaves as 6M oc (N^ — N)~^ where 
N^ — N <^ N^ is the number of orbits remaining to complete 
disruption and ^ ^ ^/(^ — !)• 

According to our estimates, tidal stripping operates for days or 
weeks, creating a relatively long-lived source of gas. The gas can 
accrete onto the black hole and produce significant electromagnetic 
radiation together with the gravitational waves observed by LISA. 
The radiation source may become bright before the mass loss spoils 
the standard EMRI template for the gravitational- wave signal. 

The gas produced by tidal stripping moves on nearly Keple- 
rian orbits that are initially close to the WD orbit. The gas probably 
leaves the star with relative velocity ^ ^;esc = {2GM/RY^'^ and 
its orbital energy differs from that of the star by a small fraction 



,/v - (Mro/MBuRy^^ - 7"'/'(M/Mbh) 
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After ^ 30 orbital periods, the differential rotation of the gas 
has stretched it into a ring around the black hole. A non-zero 
eccentricity of the donor orbit will create an eccentric ring. It 
should viscously spread and accrete onto the black hole. A small 
mass-loss fraction x can be a huge source of gas for accretion. 
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If most of the stripped matter is accreted by the black hole, the 
accretion rate is M ^ SM/P. It exceeds the Eddington value 
MEdd -- 1O^^(Mbh/1O^M0) g s"^ after N ^ 10 orbits since 
the beginning of tidal stripping, well before the final disruption of 
the WD. The accretion timescale in the viscous ring can be esti- 
mated as tacc "^ cXv'^ {H/r)~'^P where P is the WD orbital period, 
a^ = 0.01 — 0.1 is the viscosity parameter and H is the thickness 
of the ring (e.g. Shakura & Sunyaev 1973). H/r ^^ 1 is expected 
for super-Eddington accretion, which leads to tacc ^ (10 — 100) P. 

This suggests that a bright source with Eddington luminosity 
^Edd ^ 10^^ erg s~^ is created quickly, in less than 1 day after the 
beginning of tidal stripping. For a typical distance to such LISA 
sources, d ^^ 100 Mpc, the accretion ring should be detectable 
with optical and X-ray telescopes provided its approximate location 
on the sky is known. LISA is expected to localize EMRIs within 
^10 deg^ (Barack & Cutler 2004). For massive black-hole merg- 
ers, the localization information will be available weeks to months 
prior to the final coalescence (Kocsis et al. 2007, 2008; Lang & 
Hughes 2008), and the localization expectations for WD EMRIs 
are similar (S. Drasco, J. Gair, I. Mandel, E. Porter, private com- 
munications), giving sufficient time for simultaneous optical and 
X-ray observations during inspiral. 

Mass transfer in WD inspirals is special as it creates a long- 
lived source of gas very close to the black hole horizon, ro — rg ^ 
Vg. As 3. result, the donor orbit is generally not confined to a plane, 
if the black hole rotation is significant. Besides, the orbit will ex- 
perience fast precession. Therefore, the freshly stripped gas may 
collide with the previously released gas and generate shocks (e.g. 
Evans & Kochanek 1989). The resulting pattern of accretion may 
be complicated and needs careful study. 

An intriguing feature of tidal stripping is the periodic supply 
of gas. It may leave a fingerprint on the observed luminosity, mod- 
ulating it with the WD orbital period P. The modulation of M may 
create a detectable oscillation in the luminosity from the accreting 
ring, even though the accretion timescale tacc ^ (10 — 100) P ^ 
P, tends to reduce the amplitude of modulation. The shock emis- 
sion from collisions between the periodic flow from the Li point 
and the gas accumulated around the black hole may be strongly 
modulated. 

The description of mass transfer in this paper is greatly simpli- 
fied. The possibility of many repeated mass-transfer episodes is ro- 
bust, but the exact rate of tidal stripping and the dynamics of accre- 
tion need to be explored with dedicated numerical simulations. The 
great potential that such events hold for joint detections of gravi- 
tational and electromagnetic radiation provides motivation for the 
effort. The joint detection would let us witness, in real-time and 
unprecedented detail, the slow tidal stripping of a WD followed by 
its complete disruption. Note that LISA observations are expected 
to provide the mass and spin of the black hole, as well as the de- 
tails of the inspiral orbit. This can be used to model in detail the 
hydrodynamics of accretion. 

The scenario discussed in this paper assumes ro > rmin and 
uses a semi-Newtonian description for the WD orbit. A fully rela- 
tivistic model will be needed to accurately evaluate the parameter 
space for such events. The relativistic effects are especially impor- 
tant if the black hole is rapidly rotating - then ro and rmin will 
depend on the black hole mass Mbh, its spin parameter as, the 
eccentricity of the orbit, and the angle between the orbital angu- 
lar momentum and the angular momentum of the black hole. The 
competition between tidal disruption and gravitational capture by 
rotating black holes was investigated for parabolic orbits in Be- 
loborodov et al. (1992). In a broad range of Mbh, the fate of a star 



approaching the black hole depends on the orbit orientation and can 
be either disruption or capture. For inspirals with r^ < ro < rmin, 
tidal stripping does not occur. Instead, r^ reaches rmin and the or- 
bit loses stability before any mass is lost by the WD. Then the star 
is crushed by tidal forces as it falls into the black hole. This strong 
and immediate disruption is different from the gentle stripping con- 
sidered here, suggesting rich phenomenology of WD inspirals. 
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